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Abstract: We calculate the thermal spectral function of SYM plasma with finite density 
using holographic technique. The gravity dual of the finite temperature and density is 
taken as the RN-AdS black hole. In the presence of charge, linearized vector modes of 
gravitational and electromagnetic perturbation are coupled with each other. By introduc- 
ing master variables for these modes, we solve the coupled system and calculate spectral 
function. We also calculate photo-emission rate of our gauge theory plasma from spectral 
function for light like momentum. AC, dc conductivity and their density dependence is 
also computed. 
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1. Introduction 

The gauge/gravity duality ||l|, |2|, ^] opened a new possibility to quantitative study for 
strongly interacting system. Although it is not developed enough to describe realistic 
QCD, we expect to learn some features of QCD from it based on the universality of the 
hydrodynamics: in longwavelenth limit, the details of the theory does not matter. For 
example, the shear viscosity/entropy ratio [Q, ^ rj/s is universal if we neglect the higher 
derivative terms. We also expect, due to analytic structure of the theories, that there are 
similarities of super symmetric and non-supersymmetric theories can continue to the finite 
wavelenth/frequency regime. 

The quarks and gluons are liberated at high enough temperature. However, over 
Tc < T < 2 ~ 3Tc, the experimental data shows that quarks and gluons are not free but 
are strongly interacting: the small viscosity and the presence of the coherent fiow show that 
the interactions should be very strong. Such strongness of the interaction is the motivation 



why one has to abandon perturbative QCD in such energy/temperature regime. One way 
to avoid that difficulty is to rely on hQCD for the quark gluon plasma in RHIC(Relativistic 
Heavy Ion Collider). The hydrodynamic calculations of hQCD were shown to be useful to 
discuss the transport phenomena |^, ^. It is interesting to see what happens in LHC (Large 
Hadron Collider) where the energy scale is much higher p^]. 

The finite baryon density effect is very essential ingredient to understand how the core 
of neutron star and the early universe behave. It maybe uncover some significant features 
of the evolution of our universe, galaxies and stars. At RHIC experiment, the temperature 
reached is above Tc but the density is almost zero. In near future J-PARC (Japan Proton 
Accelerator Research Complex), LHC and especially FAIR (Facility For Antiproton and 
Ion Research) which probes the regime of a few times of normal nuclear density |^5| will 
tell us much about the density effect of quarks and gluons. The holographic study for the 
system with finite density in hydrodynamic regime was made in |^ |9|, |l^, 1 1 , ^ . 



The dual gravity background for the finite density and temperature is taken to be 
Reissner-Nordstrom AdS blackhole. In the phase diagram, we know much about high 
temperature, low density regime but not low temperature high density regime. Previously 
we studied some issues ie. meson mass shifted by density effect for zero temperature, finite 
density sector [g5[ . Now we will study the effect of both finite temperature and density. 
The bulk U(l) charge is identified with the particle number density of the boundary field 
theory. To see the finite frequency/momentum dependence of the response of the system, 
the spectral function is a good tool. It gives us ac conductivity and its trace is related to 
the photo-emission and di-lepton production rate [|l3| , |j, |l^. 

The spectral function with and without baryon density in the probe brane approach 



was already calculated |16|, 17, 18, 19, 21, 22, 23|. However in that approach the gravity 



back reaction to the presence of the charge is neglected. In this paper, we take bottom 
up approach where the back reaction is taked into account. We will compute the spectral 
function of tensor and vector modes which describe the fluctuation of energy momentum 
tensor and currents of hot plasma. After that the finite temperature and density effects of 
photo-emission rate are calculated and discussed. 

2. A recipe for Greens function 

In this section, we will briefiy review how to calculate thermal spectral function. To 
describe thermalized plasma holographically we need the black hole background. The 
general equations of motions for the linearized fiuctuations in this background are 

S^(n)+P(t»,q,n)^;(u) + Q(t»,q,n)^„(u)=0 (2.1) 

where E^ denote the fluctuating fields in given background and a runs 1 to n, the num- 
ber of independent fields and t», q are dimensionless frequency and momentum. Near the 
boundary, (n ~ 0), there are two local Frobenius solutions <I>i, <I>2 

^2 = u^+{l + ---) (2.2) 



01 = {i-uY 



,2 = {i-uy^f-{i + ...). 



A± is the solution of indicial equation near the boundary, A-|_ > A_ where A-j. is the 
conformal dimension of an operator and A_ is the dimension of the dual source field. Near 
the horizon, u=l, there are also two local solutions 

)i«'/2(i + ...). (2.3) 

The two different local solutions of eq.( |2.lD should be matched: 

Ea = A{w, q)$i + B{vo, q)$2 = C(rx>, q)^i + V{w, q)02. (2.4) 

However not all solutions are allowed physically because this system contains the blackhole: 
no outgoing-wave can propagate from the horizon, therefore we should impose , 2? = 0, 
which is so called the infalling boundary condition. Taking the normalization C = 1 using 
the linearity of differential equation, we have 



Mu) =^(t»,q)$i(u)+S(tx),q)<I>2(n) 



(2.5) 



Note that the coefficient A is the source of the boundary theory operator ^(t»,q) = 
Ja(n>,q), so by differentiating the generating functional twice with respect to Ja we get 
the retarded Green function. And another coefficient B corresponds to the condensate or 
vacuum expectation value of the operator Oa which couples to the source Ja- The retarded 
Green function is given by the ratio between A, B. We will give a sketch of the on-shell 
quadratic action, with eq. (|2.2| ), ( ^ ) 
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(2.6) 



obviously the first term in last line is divergent (A_ < A-|-) so it should be renormalized 
holographically p^ or we can ignore it because the imaginary part of the Greens function 
do not care about the real number which comes the first term. For the issues how to 
regulate the on-shell action, see the appendix O. The spectral function is its imaginary 



part |13|. Notice that 
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Here $« is real because the equation of motion and initial conditions are real for <l>j. Above 
expression is independent of u since it is a kind of the conserved flux P] 
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where Q = AA^, with J\f the normaUzation constant. The retarded Green function is 
defined by the recipe [0] 



= -gij{K,u = 0) i / j 



(2.9) 



From the eq. ( |2.7| ), the only thing we should know is the value of (f>i{u), ^2{u)- The spectral 
function, for example, is given when A_=0 



Xv 
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(2.10) 
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Given the normalization constant M we can calculate the spectral function numerically. 



3. RN AdS 

The dual geometry for the finite temperature and density is chosen as charged AdS black 
hole 1 12]. The action is 

S=^J d'x./^in - 2A) + ^ y d'x,/^ ^' + ^ / ^'^y^^' (3-1) 

where the cosmological constant is A = —^2F~~^^ ^^^ ^^^^ term is the Gibbons-Hawking 
term. And K is the extrinsic curvature on the boundary, 1 is the AdS radius. The metric 
of RN AdS is 



f 



dr' 



ds' = -^{-fir)dt' + dx') + ^^^^^,^ 



where the gauge charge Q is related to the black hole charge q 
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(3.2) 



(3.3) 



the five dimensional gauge theory coupling constant g^ and the gravitational constant G5 
^ can be chosen as [24| 

/ NrNf P N? 

(3.4) 
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but we will not use these parameters explicitly. The metric function /(r) is rewritten as 
fir) = ^,ir'- rl){r' - rl){r^ - rl\ r~J = ^^{l + 2 cos (^ + ^^^ (3.5) 



^Usually five dimensional gravitational constant Gs is used as k^ — SttCs but here we will use Gs as k 
itself. 



where a = +, — , and (/?+ = 47r/3, (/?__ =0, ^o = 27r/3. The charge is expressed by 6 and 



-"WflV .^a.ct.n P^firt? ^ (3.) 
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the range of outer horizon is \f^l < ?"+ < ^/ml. Finally, black hole temperature is given 
as 



47r/2 TT^n^ 2ri;^27r6V^ 2 J ^^'^^ 



where a, 6 are 

g2/2 /2 

— «-i = — 

ri 2r 



« = V' ^=T^- (3-8) 



From the horizon regularity, the black hole charge q is related with the chemical potential 

3 + 2/2^ - y/Q + 12/22 _ /x G5 ,_, 

a = Tn , where // = - -. (3.9) 

Notice that there is maximum value of q, q'^ < ■^m^l'^ which corresponds to a = 2. Horizon 
radius r^ should be real, so 1+2 cos(0/3 + 47r/3) must be positive. 

It is very useful to express the frequency and momentum as a dimesionless quantities 
vo = ■u;/(27rT),q = k/{2nT). This choice is good enough to see the finite temperature 
behavior of system but not good in zero temperature limit. Alternative way is to rescale w 
and k by the black hole radius r+, that is by b: let t» = bw,Cf = bk. At the extremal limit, 
by eq. Q, 4m^P = 27q^ and 9 = ir 

/72\l/6 //2\5/6 

r+=(y '?^/^ a = 2, b=[L^ ,-V3 (3.10) 

and the chemical potential is written 

\/3 e/ 1 

If we rescale w and k with b, in the extremal limit we rescale w and k with chemical 
potential, m ~ w/fi,q ~ k/^. 

The origin of charged black hole in string theory can be understood by STU model : 
the diagonally charged STU black hole is RN AdS. The diagonal U{1) is the subgroup of 
the SU(4) R-symmetry originally but here we assume that this U{1) is a part of flavor U(l) 
group which is relevant if we assume that the bulk filling branes||2J] are embedded in our 
AdS5 space time ^. The merit of doing this is that we can have a back-reacted gravitational 
background which is a solution of glue-quark coupled system in terms of gauge theory it 
means our approach is beyond quenched approximation. 



^Please note that this is no more than a conceptual introduction of bottom up approach in string theory. 



4. Tensor mode 

The gravitational and gauge field perturbation is classified by the boundary S0(2) ro- 
tational symmetry. This classification is summarized in the appendix ^. Tensor mode 
perturbation is easy to treat because it is completely decoupled from other fields. The 
equation of motion for hxy component is 

x// I V /) rxf I J: / „.,2 u2 . 

^5f y "•" ^4^2 



K" + '-z^K' + 1172 ( ^' - ^'/ )^« = (4-1) 



where prime denotes derivative with respect to r and hy = g^^'^^hxy Introducing new 
coordinate u = r'^/r'^, 



ds 
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{-f{u)dt^ + dx^) + ^^^——du^, f{u) = {l-u){l + u-au^). (4.2) 



u 
The equation of motion is simplified 



xll J ""/ 1X1 _| ^ / ™2 z2 . 

uf y up 



K"-'--^K' + -^A^--~^-f\hi = ^. (4.3) 



where 



'-^-2^(-3-('-f)- *-2::T(-f)-('-f)- (-) 

This differential equation has two independent solutions near boundary, 

hi = A^i{u) + B^2{u) (4.5) 

where 

$i(u) = ! + ••• + /iaog(u)$2(u), $2(u) = n2(l + •••). (4.6) 

In eq. (p.5D, we choose A = <?^i(0). By dividing A both side, we get the normalized solution 
for tensor mode 

K W = ^ = ^1 W + ^^2{u). (4.7) 



(0) ^^ ' A 



The on-shell action is given in eq. (5.7) of [[Ll| , 



S[K\ = JLj I ^A^-^Ki-k,u)K'ik,u) 
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n=0 

(4.8) 

n=l 



Our normalization is such that hy (n) — )■ 1 at the boundary. '^ By taking the imaginary part 
of the Greens function and renormalizing divergent terms, the thermal spectral function is 



""5 



x(-,c).,,, = ^^^Im 2-). (4.9) 



^More properly, we should express hy{u,k) — hy{k)hy{u), where hatted variable is the value at the 
boundary or the external source of the boundary theory. 
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Figure 1: The difference of tliermal spectral function and zero temperature one of Axxy,xy/'^- Left: 
witli fixed q =0 when /i=0(thick), 0.5(dashed), l(thin), right: Axxy,xy vs q plot and w=0.1(thick), 
0.4(dashed), 0.8(thin) with fixed ^=0.5 



Here the ratio B/A is 



B 
A 



1 
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(4.10) 
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This ratio is independent of the evaluation point. As explained before, imposing the in- 
falling condition at the horizon and Dirichlete boundary condition at the UV boundary, 
we get the numerical solution for 0i(n) and $2('u), 
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Using 6 = (1 — a/2)/27rT, one can show that the zero temperature spectral function is 

(27rT)4/3 



(4.11) 
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(4.12) 



See appendix ^ for detail. Figure ^ shows the difference between normalized thermal 



spectral function at non-zero and zero temperatures for A.Xxy,xy The thick line is zero 
chemical potential case /x=0 which is the AdS Schwarzschild case of ref. [18|. The dashed 
and solid line correspond to the fi = 0.5, 1 case respectively. When the chemical potential 
fi increases, the spectral difference grows up and the position of the small peak is shifted to 
the larger w. The position of peak in the spectral difference is the pole position of retarded 
Greens function |l^. The shift of the peak is the shift of the quasi normal mode. When 
chemical potential grows in the unit of T, the pole position grows faster than T. 

The right side of the fig. ^ shows the spectral function as a function of spatial mo- 
mentum q. When tx)=0, the peak position can be identified with the inverse screening 
length. However, for the tensor mode, there is no peak for t»=0. As we know SYM has 
conformal symmetry, hence it can not have any scale. For the finite ^, there is a broad 
peak and the peak is more sharpened when chemical potential grows. This shows that in 
a dense system the thermal particle collides more often so that the particles propagates 
shorter distance. For the light like momentum, this screening is maximized. For spacelike 
momentum, the thermal fluctuation of spectral function rapidly vanishes leaving only zero 
temperature piece. 



5. Vector mode 

Vector type perturbation consist of three independent fields, hxt,hxz,^x and equation of 
motion for these modes are coupled with each other. In hydrodynamic limit, this mode 
has a diffusion pole so that it is also named as diffusive mode. Here we are interested in 
general energy/momentum regime. The equations of motion for vector modes are 

= hf - -hf - ^^~}^ (roq/if + q^/iH - 3an5' 
= q//if + whf - 3avouB 

= hf + ^J/^hf + ^m^(u>^/.f + ro,ht) 
j/u uf^ 

= B" + Lb' + i^-l^ {vo' - qV) B-^. (5.1) 

This equation is simplified by introducing gauge invariant combination Zi = voh^ + ({hf, 
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-D + y-72 tjttA (5-2) 
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where n» = (l — |) t». It is not easy to solve these 2nd order coupled differential equations, 
but the authors of [0] have decoupled these equations by introducing master variables. 
Let us define first ^± as 

*± = -^J^Z[ + (-3au^^^^ + C±) B. (5.3) 

The equation of motion is rewritten as 



where C± is 



*i + J*i + ^^^---f l*.=0. ,5.4, 



C± = (1 + a) ± ^(1 + a)2 + 3a62A;2 = (1 + a) ± ^ (1 + a^ + 3a (l - ^ j q2 



(l+a)(l±7), where ^ = Jl + ^ ("^VqS. (5.5) 



In order to use our recipe for spectral function, we need the on-shell action for vector 
modes: 

-^-^Ki{-k,u)hl'{k,u) - 3af{u)B{-k,u)B'{k,u)) 

/ 7T^ ( ~~ ^27^1^1 - 3afBB' + — z2-:^auBht ) . (5.6) 

J (27r)4 Vuro -q/ t»-q/ / 



32Glb^ 



From the equations for master field we can get the spectral function of master fields. We 
however need the spectral function of original variables not the master fields itself. In ref. 
1 28], authors showed a systemetic way to compute the spectral function of original variables 
in terms of master variables. Let us first find the series solution of Zi,B, 

Zi = Zi(l + (a)2 - q2)^ + ...) + ^u^ + ... 
B = B{l + ---) + 7rBU + ---, 

which defines the conjugate momentums ttz^t^b- The boundary action can be written in 
terms of the boundary values of original variables and their conjugate momentums: 

l^ f d^k ( Zi-Kz 

Jbd 



32Gi64 



/(l;^(J^-'"^"^ + ---)' ^'-^^ 



where Zi , B are the boundary values of the fields and ttz,t^b are their conjugate momentum 
which will be identified with one point function and dots denote contact terms which does 
not have any derivatives with respect to u. Note that these conjugate momentums depends 
on the boundary source terms implictly. The master variables ^±1 and ^±2 have series 
solutions near the boundary, 

*± = ^±(i + --- + n±M + ---) 

Define the transformation matrix R, 

R= ( "*f^^' I, (5.8) 

the boundary value of the master fields is simply related to the boundary value of the 
original fields by R, 

-4z) <-> 





Then the conjugate momentum ttztT^b are written as 

rj\ =R-iDiag(*+,*_) r^\. (5.10) 

The boundary action is now written only by boundary values {Zi,B) and conjugate mo- 
mentum of master field II-i-. The two point function for hf and h^ is related to the gauge 
invariant variable Zi as 

_ S^Sm _ (sZiY 6^Sm _ ~2J^Sm_ ^ ( ^bY 5^Sm 1 5^Sm 



6h^6h^ yhf J 5Zi5Zi 5Zi6Zi \5A^ j 6B6B fi"^ 5B5B 

(5.11) 
therefore the correlation functions for each components are 



l^ ,nC_fl+-C+n_ 



vo 



2 



yxt,xt — Qo/^2r,4^ r^ r^ ' yxz,xz — 2 yxt,xt 

_ _,2 ^'^^ ri+-n_ 
yxt,x — yx.xt — 9 



32^5663 C+ - C- 

i c+n+-c_n_ / i/n+-n_ . . \ , ^ 

Note that when spatial momentum q or density "a" vanishes two point function Qxt,x 
vanishes. It means that the holographic operator mixing between Zi and A^ comes from the 
density effects. By following the standard recipe described in section 2, n±, the conjugate 
momentum of master fields, are computed as the ratio of two connection coefficients 



(-■) 



m±=A±[l + --- +B±u[l + --- =^± 



l + --- + Ii±u- 



(5.13) 



By comparing eq. ( 5.13| ) with eq. ( |5.13| ) we get the conjugate momentum of the master 



fields as a ratio of connection coefficient of near boundary solutions of them, 

Im -i = Im n± (5.14) 

By imposing infalling IR boundary condition for ^±, the spectral functions are computed. 

Xxtxt = 2 Im Gxtxt, Xxzxz = -^Xxtxt, Xx,xt = Im Gxxt, Xxx = 2 Im Q^x. (5.15) 

The spectral function is plotted in terms of vo, q. The figure |2| shows that the imaginary 
part of the Gxx divided by w, which is AC conductivity of thermalized plasma ( with nor- 
malization constant, 2"T2'7rT). The peak position becomes larger as the charge increases. 



The strength of that peak also increases, when charge grows. In ref.|18], they calculate only 
zero density case which is denoted by the thick line in fig. ^. The right figure in fig. y is the 
density dependence of DC conductivity. From the spectral function, DC conductivity can 
be computed by taking zero frequency limit, aE = lim^_s.o '^"''^ ~ . As density increases. 
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Figure 2: /S.Xxx/^: Deviation of finite temperature thermal spectral function from the zero tem- 
perature one, with /i=0(thick), l(dashed), 2. The normalization unit is {2'itTY^. Right : The 
density dependence of DC conductivity with normalization constant i-^27rT 



it decreases and in sufficiently large density regime, DC conductivity is negligible. This is 
rather surprizing since Drude formula in Maxwell theory says the conductivity is propor- 
tional to the density of the charge carrier. It seems that interaction between the charge 
carriers dominates the abundancy effect. Such drastic reduction of the DC conductivity 
can be an another explanation of the jet quenching phenomena which are different from 
the explanation in ref. P9| , |30|| . In highly dense system, the strongly interacting plasma 
can not carry charge over long distance because of density effect. If this is the relevant 
mechanism, raising the temperature supresses the Jet quecnching in LHC since it reduces 
Ji ~ ^l/T. 

In fig. |5|, we plot the spectral function Xxx in terms of spatial momentum with fixed 
frequency. The left one shows thick, dashed, thin line corresponds to to = 0.1, 0.3, 0.5 with 
/i=l and the right of fig. ^ also shows thick, dashed, thin line corresponds to /2=1,2,3 with 
t»=0.1. These results can be interpreted as a inverse thermal screening length of the super 
Yang-Mills plasma. It is interesting that for the tensor and vector mode the peak position 
is different. For the t» = 0, Xxx is zero so screening mass is zero. But for the finite xo there 
is the peak and the position is a function of both vo and Jl. Because the diffusive nature 
affects the interactions inside the medium, the screening effect are more complicated. 

The hydrodynamic pole in Gxx is appeared in fig at 



1 



t» ~ 



a/2 2 



2(1 + a 



(5.16) 



The left figure shows that the hydrodynamic pole position is shifted from 0.0225 (q=0.3) 
to 0.01 (q=0.2). This comes from the density effect, when ^ goes to zero the hydrodynamic 
pole in Gxx disappears, see appendix ^ This is the operator mixing result. The diffusion 
pole is only appeared in Gxt,xt or Gxz,xz not Gx,x- The right figure shows that Xxx reaches 
very rapidly to the zero temperature spectral function. 

Fig. shows real part of ac conductivity, Recr(tx)) = '^""'^ ,J. with the normalization 



unit (27rr)^Tr^. By definition, Ima 

3s 



ReG:, 



, Rea 



ImG:,^ 



The ac registivity is defined 
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Figure 3: /S.Xxx for various density, one witli a=0.5, q=0.3 (tliick;), a=0.5 q=0.2(dasiied). The 
normalization unit is {2itTY^. Right : Xj,2,(a=0.5, ct=0.3) plotted in range m £ [0,1]. 




Figure 4: Real part of AC conductivity of SYM plasma, Xxx{^,<\ = 0)/t» and the normalization 
unit is {2'kT)'^y^. Each line shows the result when p. = (thick), 1, 2, 3, 4 (thin). Left : AC 
conductivity and Right : AC registivity. 



as p{vo) = l/cr(i»). For large vo, the system has zero registivity, it means that at any density 
charge carrying is almost perfect in high frequency. 

The (xt,xt) component of spectral function Xxtxt has the diffusion pole at n? = C|^/2 
116]. The dispersion relation for diffusive channel w = Dk'^/2 gives us the diffusion constant 
D = 1/2ttT from the hydrodynamic analysis. The left figure in fig. ^ is with q=0.3 for 
various values of /2: 0.5(thick), l(dashed), 1.5(thin). When the chemical potential grows , 
the strength of peak also grows but the position itself does not. The reason of this increase 
comes from the factor q„q/2\4 i^i front of the Greens function. When p, increases the 
parameter a goes up so the overall factor (1 — a/2)~^ increases rapidly. When the system 
reaches the extremal limit, a = 2, that factor diverges and our analysis is broken down. It 
should be computed separately for the zero temperature or for the extremal RN spectral 
function from finite temperature or non extremal RN AdS blackhole. 

In the right figure of the fig. g, the peak position is shifted when q is moved. Again, 
the position is at tx) = q^/2. 
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Figure 5: The normalized thermal spectral fmiction Xxx{<\)- Left: with jl=l and varying 
t»=0.1 (thick), 0.3(dashed), 0.5(thin), Right : t»==0.1 and varying ^=0(lowest thin hne), l(thick), 
2 (dashed), 3 (solid), with normalization unit {2'kT)'^ ^^ . 
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Figure 6: The normalized thermal spectral function Xxtxt- Left: with q=0.3 and varying 
/x=0.5(thick), l(dashed), 1.5(solid), Right: with Ji—X and varying q=0.3(thick), 0.5(dashed), 
0.7(solid) with normalization unit (27rr)"^ A^„-x . 



6. Photo-emission rate 

In the heavy ion collision, the emitted photons are a good measure to see the medium 
effect. The photo emission rate of SYM plasma was calculated holographically for AdS 
Schwarzschild [13|, for D4/D8/-D8 with finite chemical potential [14| and for D3/D7 with 
finite baryon density |2^. We will focus on the photo-emission rate for our gauge theory 
dual to the RN AdS background here. Let F^ be the number of photons emitted per unit 
volume. To leading order in electromagnetic coupling e, 



dr. 



d^k 



(2vr)3 2|fc| 



^<'^C<(i^)U=fc 



C<{K) = nB(w)x^,{K) 



(6.1) 



where C^^{K) is the Fourier transformed Wightman function which is related to the spec- 
tral function multiplied by Bose-Einstein distribution function nB{w). Convert the differ- 
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(xn. 




Figure 7: The normalized trace of ttiermal spectral function Xu/"^ with light like momenta. Left 
: varying /2=0(thick), 5(dashed), lO(thin). Right : Jl— 5(dashed), lO(thin) and the plotting range 
is from zero to 30 in to. 



ential photo-emission rate into the emission rate per unit volume as a function of cj, 



~dk 



aEM 



IT 



krj^''C<{K)\^=k. 



(6.2) 



In order to calculate photo-emission rate we need to know the longitudinal part of the 
spectral function. But for the case of light like momentum uj = k, the trace of spectral 
function X/1 = v'^'^X/ii/ is obtained only by H^. From the appendix^, 



^/.^.^ret ^ (^ _ 2)n^ + n^ = 2n^ + n^ 



(6.3) 



where d is the dimension of the boundary field theory. For the light like momentum, 
the longitudinal correlator should be vanished because the projection operator diverges. 
The trace of spectral function is only given by H {w, \k\ = w). The frequency dependent 
spectral measure r|^^x^lv/^ for the light like momenta is in fig. ^. The left one in fig. ^ 
has three lines p. = O(thick), 5(dashed), lO(thin). 

In figure 7, the photo-emission rate is compared with [^], the peak position is y^max = 
1.48479/(27r) = 0.2363 the maximum value is 0.01567 with unit aEM{N^ - l)T^- The 
left figure shows the photo emission rate for /i=0(thick), l(solid), 5(dashed), lO(thin) and 
the right shows the maximum value -^{w = 0.2363) as a function of /i. Notice that 
the maximum value of the photo emission rate decreases until xx>c = 2.014 which is the 
turning point. After passing /Xc it increases. It means the thermal photon production rate 
is suppressed in low chemical potential regime but in high density regime it is enhanced. 
It may reflects the fact that at high enough density thermal screening is enhanced. As 
we have seen the suppression of conductivity and increase of thermal screening mass in 
fig. |2[ ^. When the U(l) chemical potential is very large, a ~ 2, the maximum value 
suddenly increase. Until a=l, the maximum value is slowly decreasing function of a, but 
after passing a=l it increases stiffly. But this sudden increase may be artificial, because if 
we express as the chemical potential this stiffness is not there. 
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Figure 8: The photo emission rate of SYM-EM plasma with normalization unit aEM{N^ — l)T^- 
dT^/dk with light like momenta. Left: varying /2=0(thick), l(solid), 5(dashed), lO(thin), Right : 
The maximum value of the photo emission rate as a function of chemical potential /2. 




Figure 9: The photo emission rate of hydrodynamic approximated spectral function with nor- 
malization unit {2t:T)^y^. dT^/dk with light like momenta. Left: varying /l=0(thick), l(solid), 
5(dashed), lO(thin), Right : t»^/(e^'^"^ — 1) is plotted as a function of t». 



Note that the spectral measure has maximum at the point due to the hydrodynamic 
pole of Gxx ■ 



(jt.t. 



where D 



(6.4) 



w + iDk'^' 2(1 + a) 

A simple analytic function, see appendix ^, which has only above hydrodynamic pole 
structure present the very essence of photoemission rate of dense supersymmetric Yang- 
Mills plasma. 



hydro 
A.XX 
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8eH^ 

21 
8^252 



3a 



Dk^ 



w 



1 + a D'^k'^ + ^2 



2(1 - a/2f 
(l + a)2 



hw 



2>a 



Da^^m 2(1 - a/2Y , 



+ 



vo 



Dr, 



2(1 + a) 



(6.5) 



1 + aDlq^ + w'^ (l + a)2 

Fig. m shows that the photo emission rate of hydrodynamic approximated solution. The 
left one is almost same as full numerical solution but tails in large vo are different and the 
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right one we plot t»^/(e^'^'^ — 1) as a function of to. The origin of peak in photo emission rate 
conies from the statistical factor and the change of height is density effect. This density 
effect is almost described by a single function eq. (^^) ^. 

7. Conclusion 

We have solved the equation of motion for linerized gravitational and electromagnetic 
perturbations in RN AdS background to get the holographic spectral function. Due to the 
density effect, the gravitational and electromagnetic perturbations are coupled with each 
other so it is not easy to solve. By introducing master variable, however, we can decouple 
these modes which makes the problem simpler. The problem might be handled without 
decoupling along the method discussed in [^] . 

The density effects of thermal spectral function have some interesting features. The 
boundary theory of RN AdS is believed super symmetric Yang-Mills theory with finite U(l) 
charge density. The original SYM has no dimensionful parameter and does not admit any 
quantity like screening length, energy gap, diffusion constant etc. In the finite temperature 
and density state, however, the plasma have a scale given by those. The temperature 
introduce diffusion constant, conductivity and the density modify these quantities. 

One of the interesting feature is the modification of the diffusive nature of thermal 
plasma. As we have seen the fig. § the DC conductivity is decreasing in large /i. It is 
quite natural because in dense medium, particles collide very frequently so charge carrying 
process should be suppressed. 

The photo emission rate is very important tool to probe the effects of thermal medium. 
Because photon does not interacted with other particles via strong interaction, it carries 
informations of the early stage of collision. Holographic photo emission rate is greatly 
enhanced when /2 is very large. See figure ^. Note that information of spectral function 
over only small window of w is sufficient to describe the photo emission rate, because almost 
all of the contribution for photoemission rate comes from the statistical factor l/(e'^'"^ — 1). 

It would be interesting to see whether the U(l) axial anomaly [p^, 5^, B3l can affect 
the photoemission rate. We will report this issue in the near future. 

Acknowledgments This work was supported by the WCU project of Korean Ministry 
of Education, Science and Technology (R33-2008-000-10087- 0). This work is also sup- 
ported by the National Research Foundation of Korea(NRF) grant funded by the Korea 
government (MEST) (No.2005-0049409). Also we give thanks to Yang Zhou, Alberto and 
Antonio to inform to correct some typos. 

A. The index structure of correlation functions 

The correlation functions of various operators has Lorentz index structure, and that has 
to satisfy some constraints, ie. CPT invariance. Ward identity [^]. We will briefly review 
the index structure of our correlation functions for later convenience. 

*But actually at large tt), there are deviations between fig.y andH. Please note that it is an approximated 
solution. 
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The definition of retarded correlation function of conserved current and energy mo- 
mentum tensor are 

C^,u{x -y) = -ie{x'> - y') ilMx), My)]) 
G^,a^{x -y) = -ie{x^ - y°) ([r^,(x),T,^(y)]) , (A.l) 

the state are translational invariant so we can Fourier transform these correlation functions 
into momentum space. In the equilibrium, CPT invariance told us that 

in addition correlation functions of energy momentum tensor has the property inherited 
from the symmetry of energy momentum tensor 

The Ward identity, 

k^'C^, = = k'^G^.^p (A.4) 

and if the theory has scale invariance T^ = 0, 

V^'G^.ap = 0. (A.5) 

From the Ward identity, the correlation functions are projected onto transverse spacetime 
of /c^ 

G^uo^p = P^uPo.pGb{K^) + H^,^pGs{K^) (A.6) 



where 



^^v — V^v 



K^Ky 



K2 

Hf^ual3 = -::{PiJ.aPyl3 + Pfi/sPua) — „ _ ^ PfivPaji (A- 7) 

The field theory propagator is classified by boundary S0(2) rotation symmetry |16|. 
The propagator is decomposed as scalar, vector, tensor parts according to their trans- 
formation properties under SO (2). We assume that the wave is going along z direction 
K = (t(;, 0, 0, /c) and the boundary coordinate is labeled by (t,x,y,z). The stress-energy 
tensor correlator G^yp^j is also decomposed into these categories. The conserved current, 
for our case R-current, is projected as transverse and longitudinal parts 

c^, = pj,n^ + p^^,n^ (A.8) 

where Pay,Pay is transverse and longitudinal projector which are mutually orthogonal 



flV^ flU 



T^ Ay^ f^ j nr^ n~' j nr^ , ^ , 

Pij — "ij fo~' -^mO = Pqii = 0) P^lu = ^Mi' ~ Pfiv \^-^) 
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Each component of the current-current correlator is 



C^K) = Cyy{K) = n^(i^) 



Ctt{K) 



e 



;^\ Ctz{K) 



■"^^^ cuK)=^^ii' 



For the stress-energy correlation function, the classification is slightly complicated, 

GfiualiiK) = ( P^yPaP + l^iPiiuPaP + PiiuPa/^) j ^T 

+ ( PfiuPad + 2 (-^Mi'-^a/^ + PfiiyPa/s) j ^L 



(A.IO) 



where 



^/ii^ap 



^fiuajS 



+S^yaj3Gl + Q^ual3G2 + -^/^i/a/sGs, 



~ 1 PfiaPuP + PixaPi//3 + Pfi/lPua + Pfi/sPua 



1 



L>-1 



[D 2)P^j^P^p ~'~ n — 2^>^^^'^l^ yPtiyPaP + PpiyPap, 



The transverse component of G^uaii are 

1 A;2 



G^r 



txtx 



Gx 



2 ii;2 - A;2 
1 w;2 



Gl, G; 



1 ti;/c 



txzx 



2 u;2 _ fc2 



Gi, 



1 



and longitudinal component of G are 

2G2 + 3Gl 



L^l) '^xyxy — o ^ 



Gtttt 
Gttxx 



1 



A;4 



3 (ti;2 - A;2)2 
1 fc2 



G 



1(7 



A;3 



ttt^ 



3 (ti;2 - A;2)2 



2G2 + 3Gl 



2G2 — 3G/, — 3G7^ 



(A.ll) 



(A.12) 



(A.13) 



(A.14) 



6 tt;2 — /i:2 
Note that if the theory has the scale invariance, Cl,Ct vanish. 

B. Large w spectral funtions 

In this section, we will calculate the spectral function analytically in the low and large 
frequency limit [15, ^. By using WKB methods, we can get the large frequency spectral 
function. This procedure is easily described by an example of simple Schrodinger equation, 



y"{x) + V{x)y{x) = 
here using this ansatz y{x) = e*'^^^^ we get the equation for (j){x 

-(6'f + i6" + V = 0. 



(B.l) 



(B.2) 
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Assume that (j)" is subleading, 



±\/F, 



v 



« 1^1 



(B.3) 



so first order solution is 



= ± / ^JVdx. (B.4) 

From eq. ([B.2[) we will obtain second order solution by substitute the first order solution, 



a'n2 



v±- 



i v 



±Vv + 



2VV 
iV 



4 V 



± / y/Vdx + -\nV 



so WKB solution is 



y^^) = J:^J^^SVvd.^^-^JVv<l. 



+ e" 



(B.5) 
(B.6) 



B.l Tensor mode 

For the tensor mode, the equation of motion for hy is transformed into Schrodinger form 
by choosing hy{u) = X{u)'ip{u) where X(u) = ^Ju/ f{u). 



i;" + V{u)i^ = Q, V{u) = -—, + -/-,+ 



If" , /' f" , (1-f 



7rT + 



-(tx,2-q2/) (B.7) 



4^2 4/2 2uf 2/ uf 

This equation has two singular point at u=0 and u=l. From the WKB analysis, we get 
the two linearly independent solution away from these singularities 



-01 



1 



1 



where 



VW 



p{u) 



1 - - 



) + (pi), ^2 ~ ^—-sm{S{u) + (j)2) 

Vpw 


(B.8) 


/tx)2 - q2 j^ s{u) = / p{z)dz 
Jo 


(B.9) 



Near the boundary 

Near u=0, the potential has the form 



<+|-4^ + ?l^' = » 



(B.IO) 



where Q^ = [l — ^) (ix)2 — q^). The general solution is 



^P = CiV^J2{VW^) + C2V^Y2{VW^)- 

Bessel functions have the following asymptotic forms for large x (x >> \a'^ — 1/4|) 



Ja{x) 




2 / an vr 

~ \ I — cos X 

vrx V 2 4 



), Y^{x) 




avr vr 

~ 1 / — sm I x 

vrx V 2 4 



(B.ll) 



(B.12) 



19 



We are considering large t», so asymptotic forms of Bessel is valid for our case, 

1 1 



VuJ2{V4:Q^u) ~ ^/u■ 



^cos I 2Qy/u ] 



V^y2(\/4Q2u) - v^. 



1 1 






zsm 



(B.13) 



and these are well matched with our WKB solutions p(u) ~ Q/^/u, S(u) ~ 2Q^fu. 

Near the horizon 

Near u=l, the equation of motion is 

and the solution is 



1 l + t»^ 

(l-u)2 4 



(B.14) 



^ = ^3(1 - n)i/2(i-*n') + ^4(1 - n)V2(i+*«>)_ 



(B.15) 



Near the horizon, infalling wave is only physically relevant and this boundary condition is 
to choose C4 as zero. We know the two asymptotic solutions and it should be matched at 
some point, 



1 



cos(^W-4^)+^sin(5(n)- — 



C(\-u) 



l/2-iVO/2 



from this, we get the solution for original field hy = Xip 

iJ2{^/Wu) + uY^i^fW^) = C{1 - ur''"'^ 



(B.16) 



I u.. 



(B.17) 



This is the solution for large U), q with infalling boundary condition at the black hole horizon. 
We did not fix the coefficient C yet, it is determined by the condition /i^(n = 0) = 1. Near 
the boundary, Bessel functions has series solutions and 

uJ2{^/Wu)^\q'^u + 0{u'^), uY2{^/AQ^u) L + 0(n) (B.18) 



Nothing is left, the solution is 



hi = -nQ^i uJ2{V4Qhl) + uY2{VlQh 



(B.19) 



then the spectral function is 



A-xyxy 



2ttT 



16Gl yi-a/2 



Im 






IQGl 



{w' - kye{w' - k') (B.20) 



where the theta function comes from the fact yQ^ should be real. 
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B.2 Vector mode 

For the vector mode, we have two master variables ^+, 'i>- and their equations of motions. 
These can be transformed into Schrodinger form by 

^± = -^V± (B.21) 

then equations of motion is 

The equation of motion for the master fields near the boundary is 



fe" + F±fe=0, Vi = =^L-:i;Li(l-^)- + iij-4(Ci-i- + 4^). (B^22) 



u u 



i) ~ Ci^Ji{2^/Qhi) + C2V^Yi{2^/Qhi) (B.23) 

And near the horizon, 

^ ~ C3(l - n)i/2-*«'/2 + ^(1 _ ^)i/2+*n)/2_ (g 24) 

For large Q^, two WKB solutions are 

V'l ~ J^- cos(5(tx) + </.i), ^2 ~ J-- sin(5(«) + (/-a) (B.25) 

where 

p{u) = —^Vw^-^f, S{u) = / piz)dz. (B.26) 

Vuf Jo 

This WKB solutions should be matched near boundary solutions 

2V^Jii2^/Qhl) ~ V^J ^cos (2y^Qhl-^] 

V vta/Q^u V 4 y 



2V^yi(27Q2^) ~ V^J^=sm Uy/Qhi - ^\ (B.27) 

and physically relevant boundary condition at horizon, infalling condition C4=0, 

•^ = j= = 2^c(i Ji{2^/Qhi) + Yi{2^/Qhiyj = -vrQV^A Ji{2^/Q\i) + Yi{2^^ 



(B.28) 

and by the normalization condition, ^(u = 0) = 1 the coefficient C = -ttQ/2. Then the 
spectral function for vector mode is given as 



4e262 c+ - C„ 

' .TrQ^ = J-{27TTfTT{to^-q^)e{w^-k^) (B.29) 



4g252 ^ 4g2 

Note that with Chern-Simons term, there is no difference without CS term for the large 
frequency spectral function. 
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B.3 Light like momenta 



For the light like momenta, w=k, the equation of motion is simplified |15] by 



iP±" + {w'^H + G±)7p± = 0, where H 



1-/ 

up ' 



G± 



111. I 
4/2 / 



C± 



/' , /" 



+ 



u 2 

(B.30) 

Near the black hole horizon the solution should be infalling, (1 — u)^*"2". The analytic 

solution for the equation is not known but we just consider the large to limit only to 

compute asymptotic behavior of spectral function. For large t», the leading term is tx) and 

we assume the change of the wave function in the domain u G (0, 1) is not large. 

For the large xo limit the first term in eq. (B.30) is dominant. Introduce a new variable 

c 



c 



the wavefunction is reexpressed as 



3 r 



\/—H{x)dx 



2/3 



(B.31) 



^*=li)""^ 



-H 



-1/4 



w 







(fW_ 



[m^C + l) W, where 7 = — -^ + 



/AHH"-5H'^ G 



16C^ ' V 16H» "^ ff ' ^ "'•''' 



where prime denotes derivative with respect to u, then the solution ip± is Airy function, 

^± = ^V± = ^ (^) '^' M^'/'C{n)) + ■■■= (-^) '^' Mi^y^Cin)) + ■■■ 

(B.33) 
Two point function for master variables are given as 



*'. 



± .. / 1 „ . / ?i 



II-i- = lim — ^ = lim -du In . „ , , ^.^ 

n^O^i u-,oU \l-fj Ai{vO^/^C{u)) 



C \^duAi{v,^/^C{u))^ 



(B.34) 



since / — 1 = —u'^{l + a — au), and 



lim C = ^ 9 — u + 



(B.35) 



The second term in equation is 

lim du In ( | = lim du In 

u-s>0 \1 — // «-s>0 



c 



lim du In 



u(l + a — au)y M-5>o \(l + a) 



(B.36) 



So the two point function is 
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Figure 10: The normalized trace of thermal spectral function Xa/^ with light like momenta. 
Thick line is for p, =5 and thin line is for zero temperature case. 



n± 



imn± 



Xfj, 



-,J7r/3_ 



31/3(1 + a)i/3 ( 4 + 3a + VTTa{4 + 



2/3 



{1 + VTTTly 

3^/^(1 + a)i/3 (4 + 3a + VTTa{A + 



r(3)_j^2/3 



r(^) 



2/3 



r(2/3) ^2/3 



2(i + vr+^)^ r(i/3) 

/ (27rT)2 35/6(l + a)i/3(^4 + 3a + VrT^(4 



2/3 

«) ' r(2/3) 



4e2 (1 - a/2)4/3 
C. Small to spectral funcions 



[i+vT+^y 



r(i/3) 



vo 



(B.37) 



2/3 



To check the consistency of numerical calculation, it is good to compare both large fre- 
quency and low frequency result with numerical computation. In the low frequency limit, 
authors [|ll| did the hydrodynamic analysis and it also gives us the low frequency spectral 
function. For the tensor mode, 

l^ (27rr)4 



lim Im G 



t»^.o 



xyxy 



Xxyxy 



lim 



And for the vector mode Gxx, 

_/_ 
4^ 



lim Im Gj 



3a 



l6Gl (1 - a/2)3 

/3 (27rT)^ 
IGGj (1 - o/2)3 

Dk^w 



t» 



(C.l) 



2(1 -a/2) 



(1 + a)62 D2^4 + y^2 (l + a)26 



1 . Xxx 



{2TrTf 



4e2 (1 - a/2)2 
lim lim ^^^^ = {2ttT) 



3a 



Dqk 



l + a D2q2^2 _^ „,2 



' (l + a)2 



Xxx ,_^,2 ^ 2(1 -a/2) 



t»->-oq->-o vo 



4e2 (l + a)2 ' 



(C.2) 
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where D 



2(l+a) ' 



For the Gxtxt in small w, q limit, 



ImG 



xtxt 



Xxtxt 



l^ 1 



wk 



IQGl b^ D^k^ + w;2 



/3 



{2ttTY 



t»q 



16Gi(l-a/2)3DVA:2+t»2 



(C.3) 



D. Boundary action 

In this section, we will briefly mention how to get rid of the divergences from the on-shell 
gravity action. The gauge/gravity correspondence tells us that the generating functional 
of the gauge theory is identified with the generating functional of the AdS gravity. As we 
have seen in section ||, from the generating functional we get the two point functions of the 
boundary theory G^J^J or G^n . The generating functional has some divergences which 
could be safely removed by adding counter terms, so called holographic renormalization. 
At the boundary u=0, there are two types of divergences 1/u and Log u. 

The original action eq.(p.l|) has Einstein-Hilbert, Maxwell and Gibbons-Hawking term 
and to remove the divergences we need the following counter term [p^ for the regularized 
action at the boundary. 



S. 



ct 



'ct, gravity 



s, 



ct, gauge 



1 



16GI 



d^x 



7(4) 



) + ^ loguJd^X^f^^Tmnr'^'' (D.l) 



:K] + 



where K is the curvature on the boundary. 5ct, gravity is given in ^6[. On the other 
hand. Set, gauge is obtained to cancel the logarithmic divergence coming from gauge field 
fluctuations. The boundary action for the perturbation in quadratic order derived from 



^(0) 



/3 



d^k 1 



3264^2; (2^)2^2 



uf, 
f 



^ih^r + hUh^ 



3uh' 



xl 



fK{h 



+3aBx(hf-fBj 



3uh 



xl 



(D.2) 



The Gibbons-Hawking term and the counter term ( p.l[ ) are 
l^ 1 /■ d^A; f 



5, 



(0) 



™ 32h^Glu^} (27r)2] / 



^(/i^)2-4(/i^)2-^/'(/,^)2 + 4^/,j/,f 



S. 



(0) 



3/3 



d?k 



32^4^2 n2^y (2^)2 



^\f{{K,Y-uh\{h\) 



{hir - fihir + ab'k'u'f{u)logu(B'^ 



(D.3) 

(D.4) 



24 



Then, the regularized boundary action is given as 

/3 f d'^k ( 1 



hi^^h2^{ ;('.f(-^^f'w -^j(-^^rw) 



u^0 32b^GlJ (27r)2l u 

+3aB^{-k)(hf{k) - B^'{k)) +3ab'^k'^logu{B^{-k)B^{k)) \. (D.5) 
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